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Abstract 

We report the experimental observation of the collective excitations induced in a magnetically 
trapped Bose-Einstein condensate. Low-lying mode excitations were studied by tracking the 
condensate’s center-of-mass displacement, and its aspect ratio as a function of the hold time in 
the trap. We were able to partially control the modes onset by modulating the amplitude of the 
additional field gradient used to excite the BEC. The measured excitation frequencies were found 
to be in good agreement with the literature. We have also found that the modulation amplitude 
was able to change the phase of the center-of-mass oscillation. Finally, an interesting, non-linear 
dependence was observed on the condensate aspect ratio as a function of the perturbing amplitude 
which induces the quadrupolar mode. 

PACS numbers: 03.75.Lm, 03.75.Hh 


* gugs@ifsc.usp.br 


1 




I. INTRODUCTION 


Dilute Bose-Einstein condensates (BECs) are produced and probed using atomic physics 
tool^I^, and their connection to condensed-matter physics is most clear when studying col¬ 
lective excitation^^, as well as the sound propagation^®. Much of our understanding of the 
nature of sound in quantum fluids comes from studies carried out in liquid heliumP®^. Re¬ 
cently, with gaseous condensates, the researchers have turned to these decades-old theories 
from condensed matter physics to test their validity thoroughly and very successfullj^®®. 
Nowadays, these versatile experimental system are even being used as fast complex Hamil¬ 
tonian solvers, which may provide “on demand” solutions to particular interesting problems 
in condensed matter physics. 

The nature of collective excitations in a homogeneous Bose gas depends on the hierarchy 
of three length scales: the excitation reduced wavelength, X^x', the healing length, ^ = 
l/y/STiasUo] and the mean free path for collisions, Xmfp- In a inhomogeneous Bose trapped 
gas, the excitations are influenced by an additional length: the typical condensate size, the 
Thomas-Fermi radii. This gives rise to three different excitation regimes: Hydrodynamic, 
collision and collisionless. The hrst observations of the low-lying excitations in gaseous 
BECs were made at JILA® and at the MIT® in the zero temperature limit. Soon after, 
Stringari carried on in depth theoretical studies of the collective modes creating the 

base framework to understand the observations. 

In this work, we investigate the collective excitations of ®^Rb condensates held in a cigar¬ 
shaped Quadrupole-Ioffe conhguration (QUIC) trapP^. The collective excitations of the con¬ 
densate are produced by an additional, time varying, magnetic held gradient superimposed 
to a QUIC potential. Besides the expected oscillation of the center of mass, we observed 
the BEC undergoing shape oscillations. The frequency of the observed collective excitations 
were measured and compared to the theoretical predictions. We see evidence of a non-linear 
dependence on the amplitude of the perturbing held gradient, not yet predicted by the 
theory. We believe that the study of collective modes is relevant because they present the 
fundamental features of quantum degenerate gases. Furthermore, we would like to better 
understand the vortex nucleation and the onset of turbulence in perturbed condensates. 

In the next sections, we present a brief discussion on the collective modes observed in 
BECs, followed by a description of the experimental setup used to acquire the data. The 
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results are then presented and discussed based on the current theory for the dispersion 
relations, where good agreement is found for the measured mode frequencies. Finally, the 
conclusions are presented. 

II. COLLECTIVE EXCITATIONS IN GASEOUS BECS 

Bose-Einstein condensates were object of intensive studies early since they were hrst pro¬ 
duced in gaseous alkali metal samples. The low-lying collective excitations were theoretically 
investigated and predicted just one year after the initial experimental BEC observation^^*^. 
Back then, StringarP was able to derive analytic solutions for both isotropic and anisotropic 
harmonic potentials in the Thomas-Fermi regime. 

The normal modes of a BEC are classihed by quantum numbers {n,l,m), where n is 
the radial quantum number and I, m stand for the total angular momentum and its axial 
projection quantum numbers, respectively. For cylindrical symmetry (as in our trap), m is 
still a good quantum number, but I is not. Thus, the normal modes are a superposition 
of wave functions with the same m. Stringari discussed an anisotropic harmonic potential 
with axial symmetry and showed that the lowest m = 0 modes are coupled excitations of 
the (0, 2, 0) modes, quadrupolar surface oscillation, and also the (1, 0, 0) symmetries. For a 
cigar-shaped condensates the lowest modes were predicted at frequencied^^^^® 

uP{m = ±/) = IujI, ( 1 ) 

a;^(m = ±/ — 1) = (/ — l)a;^ -f ojI, (2) 

and 

= 0) = ul{2 + T - 16A2 + 16). (3) 

In the above equations Ur are the axial, radial trap frequencies which will be properly 
dehned in the next section, and A is the ratio of those two frequencies known as the asym¬ 
metry parameter. In the large anisotropy limit (A << 1), the solutions of Eq. above are: 

and The hrst, low-lying, shape oscillations associated with the above disper¬ 
sion relations are: 1) the cylindrically symmetric (slow) m = 0 quadrupole mode where the 
axial (z—axis) length oscillates out of phase with the radial lengths; 2) the in-phase (fast), 
mostly radial, m = 1 breathing mode; and 3) the non-symmetric radial oscillation m = 2 
quadrupole mode where the two radial lengths oscillate out of phase. 
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III. EXPERIMENTAL SETUP 


The experimental sequence to produce the BEG runs as follows. First, BECs con¬ 
taining about 1 — 2 X 10^ atoms are formed in the \F = 2,mF = +2 > hyperhne state, 
with a small thermal fraction in a cigar-shaped QUIC magnetic trap. In order to achieve 
the degeneracy, the atoms undergo forced evaporation for 23s, going from 20MHz down to 
about 1.71MHz, following a non analytic curve experimentally determined by optimizing 
the phase space density and the runaway condition in each step. We typically produce sam¬ 
ples with no more than 30 or 35% of thermal fraction, well into the Thomas-Fermi regime. 
The measured trapping frequencies are: = 27r x 21.1(l)Hz in the symmetry axis, and 

Ur = 271 X 188.2(3)Hz in the radial direction. The radial trapping frequency depends on the 
offset held, Bq, at the trap bottom and the magnetic held gradient in the radial direction, 
which can be written as: 



where mp, Hb and gp are the magnetic quantum number, the Bohr magneton, and Lande 
g-factor for the hyperhne state \F = 2,mp = +2 >, and m is the atomic mass. The axial 
trapping frequency is determined by the axial curvature B” through Uz = ^/mpgpfXBB”/m. 
The values of Bq, B', and B" are determined based on the geometric parameters of the QUIC 
coils and the electric current running through them. For our system, the calculated values 
are Bq = 1.22G, B'^ = 162G/cm, R' = 165G/cm, and B” = 282.8G/cm^. Correspondingly, 
the theoretical results of the radial and axial trapping frequencies are = 27r x 188.6Hz 
and Uz = 277 X 21.4Hz, respectively, for a current of 25A through the coils. Finally, we have 
mp = 2, gp = 1/2, and the anisotropy parameter is A = ^ = 0.113. 

To induce the collective excitations, an additional sinusoidal magnetic held gradient is 
turned on and superimposed on the atoms just after the BEC is produced. The external gra¬ 
dient is a quadrupolar held generated by an anti-Helmholtz coil pair set with symmetry axis 
closely parallel to the weak trap axis. The coil pair is not perfectly aligned to the condensate 
axis {6 < 5°), so its held gradient, near the trapping region, has components parallel to each 
trap eigen-axis. A superposition of shape morphing, trap bottom displacement and axial 
rotation is coupled to the BEC. Additional, sinusoidal held gradients going from zero up to 
800mG/cm, along the vertical axis, were applied to the trapped atoms inducing the onset 
of the collective modes. The extra gradient oscillating at 189Hz was kept on for six extra 
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periods of trap oscillation and then switched-off. In Fig{T} we present a typical set of the 
perturbed BEC samples imaged for different hold times, nearly completing a full oscillation 
cycle. 

We then kept the BEC held in the trap for different time intervals, thoid, ranging from 
30ms to 50ms after which we observe it by performing a resonant time-of-flight (TOE) 
absorption imaging, after 20ms free fall to avoid pixel saturation and increase the accuracy. 
Center-of-mass position of the cloud as well as its axes sizes and orientations are recorded 
and plotted as a function of the thoid, from where we extract both frequency and amplitude 
of the collective modes. 


IV. RESULTS AND DISCUSSIONS 

Using the technique explained in the previous section, we are able to observe dipolar, 
quadrupolar, breathing and also the scissors modes in our BECs. For small excitation 
amplitudes, we observed the tilting about the gravity direction {y — axis). Such periodic 
axial tilting is generated by a sudden rotation of the conhning trap, and it is known in the 
literature as scissors mode. It was theoretically studied by Query- Odeliii^ and others and 
will not be discussed here. Here, we will focus on the dipolar and the quadrupolar modes 
only. 

After being perturbed, the BEC center-of-mass always oscillate in the dipolar mod^, 
into the trap. This mode corresponds to the full oscillation at the trap frequency in each 
direction. We have also observed the onset of the quadrupolar mod^. Since the longitudinal 
and transverse directions do not oscillate exactly at 90° out-of-phase, we believe that the 
higher order modes may also be onset. Moreover, the condensate’s long axis cyclic tilting 
was observed and priviously reported^, see Figjl} The observation of the scissors modd^ 
shows that the perturbing held gradient is able to cause the BEC to rotate. 

The study of the low-lying collective excitations was carried on by following the center- 
of-mass (CM) motion and the aspect ratio (AR) of the BECs as a function of the hold 
time. For small amplitudes, the CM is usually located near the center of the atomic density 
distribution, which is also the geometric center of the BEC. But, in some situations, at large 
perturbing amplitudes (> 80%), we observed slosh dynamics taking place, clearly displacing 
the thermal and condensed fractions in opposite, out-of-phase, CM motionl^. It is important 
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to mention that the dipolar mode preserves the AR when compared to an unperturbed BEC. 
Therefore, the AR analysis are used only to study the higher order modes (quadrupolar and 
breathing). 

In Figg a), we present the results acquired by directly tracking the center-of-mass (CM) 
motion as a function of the hold time for four different perturbing amplitudes. One may note 
that, the perturbing energy not only increases the CM displacement from the rest position, 
but also shifts the phase of the sinusoidal motion. Figgb) presents the AR evolution with 
the hold time, measured for the quadrupolar mode. Again, note the phase shift introduced 
in the AR evolution as the perturbing energy amplitude raises from 80mG/cm (10%) to 
553mG/cm (70%). The data shown in Figsga) and (b) are htted by simple sinusoidal 
functions from where we determine the oscillation frequency of the collective modes and 
their amplitudes. 

The evolution of the CM and the AR as a function of the perturbing amplitude is shown 
in Figga). On the one hand, the CM oscillation amplitude proportionally increases with 
the excitation amplitude starting from zero, as expected. On the other hand, the evolution 
of the AR clearly show different evolution slopes, depending on perturbing amplitude range. 
And it displays an interesting non-linear dependence near 400mG/cm. A small region around 
350 — 400mG/cm presents a interesting non-linear behavior, rapidly evolving from the initial 
to the hnal (steeper) slope, strongly suggesting the emergence of a different regime for 
perturbing amplitudes lager than 400mG/cm. We do not have yet a theoretical explanation 
for this observation. 

Figgb) presents the changes on the dipolar and quadrupolar oscillation frequencies as 
the perturbing amplitude increases. We measured Ud = 27r x 185Hz for the dipolar oscilla¬ 
tions, compared to Ur predicted. And, ojq = 2 ti x 375Hz = for the fast oscillation, 

quadrupolar mode, which is to be compared to the theoretical prediction: 2ur. The ob¬ 
served modes were previously studied and analytic expressions for the dispersion relations 
were theoretically derivecP®. The agreement found between our results and the predicted 
frequencies, Eqsj^ to presented in section 2 is very good. The results show almost no 
change with the frequency for amplitudes up to about 400mG/cm and a slight decrease 
(2 — 3%) for larger amplitudes. Similar observations were seen befor^. We expected to see a 
constant evolution along the full range and believe the slight change may attributed to non¬ 
linear interactions coupling the condensate normal mode^. For nearly pure condensates. 
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the thermal excitations should have no effect at all and, to the best of our knowledge, there 
is still no theoretical explanation for a possible slow down of the collective modes frequency 
of a trapped BEC. 

We believe though, that the different behavior observed in between the CM oscillation and 
the shape oscillation (quadrupolar) might be due to way we perturb the BEC. The addition 
of a small amplitude, time-varying, tilted magnetic field gradient will mostly displace the 
CM of the BEC. But, once the axial and the radial trapping gradients are a factor of two 
different, it is enough to induce rotation and shape oscillations. This observation is in 
agreement with results presented in Figsj^-|^ Therefore, small perturbation results in CM 
displacement with little effect on the BEC shape. 

As the perturbation amplitude increases, the condensate is pushed farther away from 
the trap bottom. In a QUIC trap, the BEC needs to be displaced away from the trap 
center, up to regions with steeper field curvature (see Fig.2(d) in the original papei^, 
where the unbalanced compression is likely to be effective inducing shape oscillations. That 
is consistent with the near resonant perturbation used to drive our BEC. Moreover, the QUIC 
traps are known to be asymmetric along the weak axis direction due to the presence of a 
single Ioffe coil, making it easier for the condensate to reach the non symmetric boundaries 
of the trapping potential. This also seems consistent with the observation of a change in the 
curvature shown in Figj^ where the BEC breaks free from the inner harmonic (symmetric) 
trapping region and travels through the outermost, anharmonic, trap region. These outer 
trap regions present much stepper field gradient and curvature making them suitable to 
produce shape oscillations in a BEC by unbalanced compression. Beyond that, we believe 
that the steep gradient regions might be the inducing mode coupling and and even the 
dephasing observed in Figj^ 

This complex coupling dynamics seems particular to our experimental system. It may 
also be the key to the better understand the vortex nucleation and the turbulence emergence 
observed in previous experiment^. In contrast, the presented results are in good agreement 
with the theoretical predictions except for the observed threshold for the quadrupolar shape 
oscillations. 
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V. CONCLUSIONS 


We have studied the excitations induced in cigar-shaped condensates produced in a QUIC 
trap, near zero temperature. Collective modes were excited by sinusoidally varying currents 
running through additional coils, externally set for this purpose. The excitation was not 
axially symmetric, and so the m 7^ 0 modes were excited too. Dipolar oscillations were 
readily induced and used to accurately measure the trap frequencies. At least, two additional 
shape oscillations were set on and studied. The higher frequency, primarily, radial breathing 
mode was excited and its frequency was measured, in very good agreement with the predicted 
theoretical value. Finally, we have observed an interesting, nonlinear, onset dependence on 
the minimum excitation amplitude needed to start changing aspect ratio. 
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FIG. 1. Typical dataset acquired during the experimental runs. It shows the collective excitations 
introduced by the additional magnetic field gradient, for six different hold times, before the trap 
was switched-off. Please, note the dipolar, quadrupolar, breathing and scissors modes, all coupled, 
appearing in the normalized absorption image snapshopts taken after 20ms time-of-flight. The 
perturbing amplitude increase mostly stretches the BEC radii and displaces its center of mass 
farther from the rest position (no perturbation^. 










FIG. 2. (Color online) (a) the center-of-mass (CM) motion as a function of the hold time for differ¬ 
ent perturbing amplitudes, (b) The AR amplitude versus the hold time, for 4 different amplitudes 
from 10 to 70% of 800mG/cm, corresponding to the quadrupolar collective mode. 
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FIG. 3. (Color online) (a) CM displacement and the AR evolution as a function of the perturbing 
gradient absolute values, (b) Dipolar and quadrupolar oscillation frequencies versus the perturbing 
amplitude. 
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